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O ■ Abstract 

This is a very brief introduction to Wilson's Renormalization Group with emphasis on math- 
. ematical developments. 

'Qh 1 1 Introduction 



The renormalization group (RG) in its modern form was invented by K.G. Wilson in the context of 
statistical mechanics and euclidean quantum field theory (EQFT). It offers the deepest understand- 
ing of renormalization in quantum field theory by connecting EQFT with the the theory of second 
order phase transition and associated critical phenomena. Thermodynamic functions of many sta- 
^ . tistical mechanical models ( the prototype being the Ising model in two or more dimensions) exhibit 

power like singularities as the temperature approaches a critical value. One of the major triumphs 
of the Wilson RG was the prediction of the exponents ( known as critical exponents) associated to 
these singularities. Wilson's fundamental contribution was to realize that many length scales begin 
to cooperate as one approaches criticality and that one should disentangle them and treat them one 
at a time. This leads to an iterative procedure known as the renormalization group. Singularities 
and critical exponents then arise from a limiting process. Ultraviolet singularities of field theory 
can also be understood in the same way. Wilson reviews this in [WK 1974] and gives the historical 
genesis of his ideas in [W 1983]. 

The early work in the subject was heuristic in the sense that clever but uncontrolled approx- 
imations were made to the exact equations often with much success. Subsequently authors with 
mathematical bent began to use the underlying ideas to prove theorems: Bleher and Sinai [BS 1975] 
in work related to critical phenomena and Benfatto, Cassandro, Gallavotti et al [BCGNOPS 1978, 
1980] in the construction of super-renormalizable quantum field theories. The subject saw further 
mathematical development in the work of Gawedzki and Kupianen [ GK 1980, 1983, 1984, 1985] 
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and Balaban [1982a, 1982b]. Balaban in a series papers ending in [B 1989] proved a basic result 
on the continuum limit of Wilson's lattice gauge theory. Brydges and Yau [BY 1990] simplified 
the mathematical treatment of the renormalization group for a class of models and this has lead to 
further systemization and simplification in the work of Brydges, Dimock and Hurd [BDH 1998] and 
Brydges, Mitter and Scoppola [BMS 2003]. Another method which has been intensely developed 
during the same historical period is based on phase cell expansions: Feldman, Magnen, Rivasseau 
and Seneor [FMRS 1986, 1987] developed the early phase cell ideas of Glimm and Jaffe and were 
able to prove independently many of the results cited earlier. Although these methods share many 
features of the Wilson RG they are different in methodology and thus remain outside of the purview 
of the present exposition. 

A somewhat different line of development has been the use of the RG to give simple proofs of 
perturbative renormalizability of various quantum field theories: Gallavotti and Nicolo [GN 1985] 
via iterative methods and Polchinski [P 1984] who exploited a continuous version of the RG for 
which Wilson [W 1974] had derived a non-linear differential equation. Polchinski's work has been 
extended to non-abelian gauge theories by Kopper and Muller [KM 2000]. 

Finally it should be mentioned that apart from quantum field theory and statistical mechanics 
the RG method has proved fruitful in other domains. An example is the study of interacting 
fermionic systems in condensed matter physics ( see the articles of Mastropietro and Salmhofer in 
this volume). In the rest of this article our focus will be on quantum field theory and statistical 
mechanics. 

2 The RG as a Discrete Semigroup 

We will first define a discrete version of the RG and consider its continuous version later. As we 
will see the RG is really a semigroup. Thus calling it a group is a misnomer. 

Let <j> be a Gaussian random field (see e.g. Gelfand and Vilenkin, [GV 1964] for random fields) 
in R d . Associated to it there is a positive definite function which is identified as its covariance. In 
quantum field theory one is interested in the covariance 

E{4>(x)4>(y)) = const.\x - y\- 2 ^ = dp j*<*~v) (2.1) 

Here [(f)] > is the (canonical) dimension of the field, which for the standard massless free field is 
[(f)] = The latter is positive for d > 2. However other choices are possible but in EQFT they 
are restricted by Osterwalder-Schrader positivity. This is assured if [4>] = with < a < 2. If 
a < 2 we get a generalized free field. 

Observe that the covariance is singular for x = y and this singularity is responsable for the 
ultraviolet divergences of quantum field theory. This singularity has to be initially cutoff and there 
are many ways to do this. A simple way to do this is as follows. Let u(x) be a smooth, rotationally 
invariant, positive definite function of fast decrease. Examples of such functions are legion. Observe 
that 

\ x - y\~^ = CO n S t. [°° ^ r 2 ^ u(^) (2.2) 
Jo I I 
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as can be seen by scaling in I. We define the unit ultraviolet cutoff covariance C by cutting off at 
the lower end point of the I integration ( responsable for the singularity at x = y) at I = 1, 



OO 



C(x-y)=J i fr^u^) (2.3) 

C(x — y) is positive definite and everywhere smooth. Being positive definite it qualifies as the 
covariance of a Gaussian probability measure denoted fie on a function space f2 ( which it is not 
necessary to specify any further). The covariance C being smooth implies that the sample fields of 
the measure are fie almost everywhere sufficiently differentiable. 



Remark : Note that more generally we could have cutoff the lower end point singularity in (j2.1|) at 
any e > 0. The e-cutoff covariance is related to the unit cutoff covariance by a scale transformation 
( defined below) and we will exploit this relation later. 

Let L > 1 be any real number. We define a scale transformation Sl on fields <p by 

S L <f>{x) = L-M 0(|) (2.4) 



on covariances by 



S L C(x-y)=L- 2 ^C(^-) (2.5) 



and on functions of fields F(4>) by 



S L F(<P) = F(S L <P) (2.6) 



The scale transformations form a multiplicative group : S"£ = Sl™- 
Now define a fluctuation covariance Tl 

r L {z-y) = £frMu£Z*) (2.7) 

— v) is smooth, positive definite and is of fast decrease on scale L. It generates a key scaling 
decomposition 



C(x -y)= T L (x - y) + S L C(x - y) (2i 



Iterating this we get 



where 



C(x - y) = J2 T n(x - y) (2.9) 

71=0 

T n (x -y) = S L nT L (x -y) = L~ 2n ^ T L {^-) (2.10) 

The T n (x — y) are of fast decrease on scale L n+l . 

Thus (|2,9|) achieves the decomposition into a sum over increasing length scales as desired. Being 
positive definite the T n qualify as covariances of Gaussian probability measures , and therefore 



3 



A*C = (£) Mr„- Correspondingly introduce a family of independent Gaussian random fields Q n , 

n=0 

called fluctuation fields, distributed according to /ir„ • Then 

oo 

<t> = Cn (2.H) 

n=0 

Note that the fluctuation fields £ n are slowly varying over length scales L n . In fact an easy estimate 
using a Tchebycheff inequality shows that for any 7 > 



y\ < L n ^ fi C (\(n(x) - Cn(y)\ > 7) < const.7- 2 (2.12) 



which reveals the slowly varying nature of £ n on scale L n . ()2.11|) is an example of a multiscale 
decomposition of a Gaussian random field. 

The above implies that the \iq integral of a function can be written as a multiple integral over 
the fields ( n . We calculate it by integrating out the fluctuation fields ( n step by step going from 
shorter to longer length scales. This can be accomplished by the iteration of a single transformation 
Tl, a renormalization group transformation, as follows. Let F(cj)) be a function of fields. Then we 
define a RG transformation F — > T^F by 

(T L F)(</>) = S l h Tl * F(<f>) = J durAQ F(( + S^) (2.13) 

Thus the renormalization group transformation consists of a convolution with the fluctuation mea- 
sure followed by a rescaling. 

Semigroup property : The discrete RG transformations form a semigroup 

T L T Ln = T Ln+ i (2.14) 

for all n > 0. 

To prove this we must first see how scaling commutes with convolution with a measure. We 
have the property : 

/Ur L * S L F = S Lf is L r L * F (2.15) 

To see this observe first that if £ is a Gaussian random field distributed with covariance Tl then 
the Gaussian field Sl( is distributed according to Sl^l- This can be checked by computing the 
covariance of SlC- Now the lefthand side of (|2.15j) is just the integral of F(SlC + Sl4>) with re- 
spect to dfir L (C)- By the previous observation this is the integral of F(( + Sl4>) with respect to 
dfis L r L ((), and the latter is the right hand side of (|2,15j) . 

Now we can check the semigroup property trivially : 



T L T L nF = Sl^f l * Sl^^t L u * F = SlSl^^s L uT l * Hr L n * F 
= S Ln +ifir Ln +s L nr L * F = S L n+ifJ,r rn+1 * F 
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(2.16) 



We have used the fact that T^n + Sl^l = T L n+i. This is because Sl^l has the representation 
(|2.7|) with integration interval changed to [L n , L n+l ]. 



We note some properties of Tl- Tl has an unique invariant measure, namely nc ■ F° r an Y bounded 
function F 

j d^cT L F = J dfi c F (2.17) 

To understand (|2.17|) recall the earlier observation that if <f> is distributed according to the covariance 
C, then Sl4> is distributed according to SlC. By (|2.8|) Tl + SlC = C. Therefore 



diL C T L F = j dfj, c SLl^T L * F 
dlis L cl*V L * F 

dfi c F (2.18) 



The uniqueness of the invariant measure follows from the fact that the semigroup Tl is realized by 
a convolution with a probability measure and therefore is positivity improving: F > 0, \iq a.e. 
T L F > 0, nc a.e. 

Finally note that Tl is a contraction semigroup on L p (d^c) f° r 1 < P < °°- To see this note 
that since Tl is a convolution with a probability measure : TlF = /is _ir L * SlF we have via 
Holder's inequality \TlF\ p < Tl\F\ p . Then use the fact that \ic is an invariant measure. 



Eigen functions : Let : p nt m '■ {^>{x)) be a C Wick ordered local monomial of m fields with n 
derivatives. Define P n . m {X) = f x dx : p njm '-c ( x )- The P njm (X) play the role of eigenfunctions of 
the RG transformation Tl upto a scaling of volume: 

T L Pn,m(X) = F d - m ^- n P n ^{L- l X) (2.19) 

Because of the scaling in volume the P n}Tn (X) are not true eigenfunctions. Nevertheless they are 
very useful because they play an important role in the analysis of the evolution of the dynamical 
system which we will later associate with Tl- They are classified as expanding (relevant), contract- 
ing (irrelevant) or central (marginal) depending on whether the exponent of L on the right hand 
side of 12.191 is positive, negative or zero. This depends of course on the space dimension d and the 
field dimension [(f)]. 

Gaussian measures are of limited interest. But we can create new measures by perturbing the 
Gaussian measure \ic with local interactions. We cannot study directly the situation where the in- 
teractions are in infinite volume. Instead we put them in a very large volume which will eventually 
go to infinity. We have a ratio of two length scales, one from the size of the diameter of the volume 
and the other from the ultraviolet cutoff in [ic and this ratio is enormous. The Renormalization 
group is useful whenever there are two length scales whose ratio is very large. It permits us to 
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do a scale by scale analysis and at each step the volume is reduced at the cost of changing the 
interactions. The largeness of the ratio is reflected in the large number of steps to be accomplished 
this number tending eventually to infinity. This large number of steps has to be mathematically 
controlled. 



Perturbation of the Gaussian measure 

Let Aat = [-^r,^-] d C H d be a large cube in R d . For any X C A N let Vq(X,</>) be a local 
semi bounded function where the fields are restricted to the set X. Local means that if X, Y are 
sets with disjoint interiors then Vq{X \JY,(f>) = Vq(X) + Vq(Y). Consider the integral ( known as 
the partition function in QFT and statistical mechanics) 

Z(A N ) = J dfx c ((t>)z (A N ,ct>) (2.20) 

where 

zo(X,<f>) =e~W0 (2.21) 

and 

V 0) (Aiv, 0) = d» c (<t>) e- y °CW) (2.22) 

is the corresponding probability measure. Vq is typically not quadratic in the fields and therefore 
leads to a non-gaussian perturbation. For example , 



V (X,, 



-- [ dx{t\V<i){x)\ 2 + g <j)\x)+i M) <p 2 {x)) (2.23) 

where we take go > 0. The integral (|2.2U|) is well defined because the sample fields are smooth. 

We now proceed to the scale by scale analysis mentioned earlier. Because fic is an invariant 
measure of Tl we have for the partition function Z(An) in the volume A at 

Z(A N ) = J d/ic(<f>)zo(A N ,<l,) = J ducifiTLZaiAN,® (2.24) 

The integrand on the right hand side is a new function of fields which because of the final scaling 
live in the smaller volume Ajv-i- This leads to the definition : 



zi(Ajv_i, 0) = T L z (A N , 0) (2.25) 

Because Vq was local zq has a factorization property for unions of sets with disjoint interiors. This 
is no longer the case for z\. Wilson noted that nevertheless the integral is well approximated by 
an integrand which does but the approximator has new coupling constants. The phrase "well ap- 
proximated" is what all the rigorous work is about and this was not evident in the early Wilsonian 
era. The idea is to extract out a local part and also consider the remainder. The local part leads 
to a flow of coupling constants and the (unexponentiated ) remainder is an irrelevant term. This 
operation and its mathematical control is an essential feature of RG analysis. 
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Iterating the above transformation we get for all < n < N 

z n+ i(Ajv-n-l, <j>) = T L z n (A N ^ n , 4>) (2.26) 

After N iterations we get 

Z(A N ) = I dii C {<f>)zN(Ac,4>) (2-27) 



where Ao is the unit cube. To take the N — > oo we have to control the infinite sequence of itera- 
tions. We cannot hope to control the infinite sequence at the level of the entire partition function. 
Instead one chooses representative coordinates for which the infinite sequence has a chance of hav- 
ing a meaning. The coordinates are provided by the coupling constants of the extracted local part 
and the irrelevant terms, (an approximate calculation of the flow of coupling constants is given in 
Section 3). The existence of a global trajectory for such coordinates helps us to control the limit 
for moments of the probability measure (correlation functions). The question of coordinates and 
the representation of the irrelevant terms will be taken up in Section 4. 

Ultraviolet cutoff removal : 

The next issue is ultraviolet cutoff removal in field theory. This problem can be put into the 
earlier framework as follows. Let e^v be a sequence of positive numbers which tend to as N — > oo 
. Following the Remark after (|2.3j) we replace the unit cutoff covariance C by the covariance C tN 
defined by taking ejy instead of 1 as the lower end point in the integral (|2.3|) . Thus acts as a 
short distance or ultraviolet cutoff. It is easy to see that 



Consider the partition function Z 6N (A) in a cube A = r — — — ld • 



C €N (x-y) = S eN C(x-y) (2.28) 

R R 
2 ' 2 J • 

Z eN (A) = J dii CiN {<t>) e- v °( A >*'t»,§»,w) (2.29) 
where Vq is given by (|2.23|) with go, Ho replaced by gNifiN- By dimensional analysis we can write 

t AT _ M4>]-d+2)t 3 „_,(4M-d) ,- J „_,(2M-d) #< (nor* 
t,N — % S) 9n — g, UN — e N A* (^.OUJ 

where g,£,n are dimensionless parameters. Now <f> distributed according to C eN equals in distribu- 
tion S eN 4> distributed according to C. Therefore choosing ejy = L~ N we get 



Z eN (A)= I dnc(<j>) e - y o(A,S, Ar *,eAr,gAr.A*r) 

duci^e-v^^^g,*) (2.3i) 



where Am = [—L N ^, L N ^] d . Thus the field theory problem of ultraviolet cutoff removal i.e. taking 
the limit eAr — > has been reduced to the study of a statistical mechanical model in a very large 
volume. The latter has to be analyzed via RG iterations as before. 
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Critical field theories 

As mentionned earlier we have to study the flow of local interactions as well as that of irrelevant 
terms. Together they constitute the RG trajectory and we have to prove that it exists globally. 
In general the trajectory will tend to explode after a large number of iterations due to growing 
relevant terms (characterized in terms of the expanding Wick monomials mentionned earlier). 
Wilson pointed out that the saving factor is to exploit fixed points and invariant manifolds of 
the fixed points by tuning the initial interaction so that the RG has a global trajectory. This 
leads to the notion of a critical manifold which can be defined as follows. A fixed point will 
have contracting and/or marginal attractive directions besides expanding ones. In the language of 
dynamical systems the critical manifold is the stable or center stable manifold of the fixed point 
in question. This is determined by a detailed study of the discrete flow. In the examples above it 
amounts to fixing the initial "mass" parameter = Mc(<7o) with a suitable function \i c such that 
the flow remains bounded in an invariant set. The critical manifold is then the graph of a function 
from the space of contracting and marginal variables to the space of /u's which remains invariant 
under the flow. Restricted to it the flow will now converge to a fixed point. All reference to initial 
coupling constants have disappeared. The result is known as a critical theory. 

Critical theories have been rigorously constructed in a number of cases. Take the standard (j) 4 
in d dimensions. Then [(f)] = For d > 5 the </> 4 interaction is irrelevant and the Gaussian fixed 
point is attractive with one unstable direction ( corresponding to fi ). In this case one can prove that 
the interactions converge exponentially fast to the Gaussian fixed point on the critical manifold. 
For d = 4 the interaction is marginal and the Gaussian fixed point attractive for g > 0. The critical 
theory has been constructed by Gawedzki and Kupainen in [GK 1985] starting with a sufficiently 
small coupling constant. The fixed point is Gaussian ( interactions vanish in the limit) and the 
convergence rate is logarithmic. This is thus a mean field theory with logarithmic corrections as 
expected on heuristic grounds. The mathematical construction of the critical theory in d = 3 is an 
open problem. ( It is expected to exist with a non-Gaussian fixed point and this is indicated by the 
perturbative e expansion of Wilson and Fisher in 4 — e dimensions ) . However the critical theory 
for d = 3 for [</)] = 2=£ for e > held very small has been rigorously constructed by Brydges, Mitter 
and Scoppola in [BMS 2003]. This theory has a nontrivial hyperbolic fixed point of O(e). The 
stable manifold is constructed in a small neighborhood of the fixed point. Note that the ( uncutoff ) 
covariance is Osterwalder-Schrader positive and thus this is a candidate for a non-trivial EQFT. 
For e = 1 we have the standard situation in d = 3 and this remains open as mentioned earlier. A 
very simplified picture of the above is furnished by the perturbative computation in Section 3. 

Unstable fixed points 

We may attempt to construct field theories around unstable fixed points. In this case the initial 
parameters have to be adjusted as functions of the cutoff in such a way as to stabilize the flow in the 
neighborhood of the fixed point. This may be called a genuine renormalization. A famous example 
of this is pure Yang-Mills theory in d = 4 where the Gaussian fixed point has only marginal unstable 
directions. Balaban in a series of papers ending in [B] considered Wilson's lattice cutoff version of 
Yang-Mills theory in d = 4 with initial coupling fixed by the two loop asymptotic freedom formula. 
He proved by lattice RG iterations that in the weak coupling regime the free energy per unit volume 
is bounded above and below by constants independent of the lattice spacing. Instability of the flow 
is expected to lead to mass generation for observables but this is a famous open problem. Another 
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example is the standard non-linear sigma model for d=2. Here too the flow is unstable around the 
Gaussian fixed point and we can set the initial coupling constant by the 2-loop asymptotic freedom 
formula. Although much is known via approximation methods ( as well as by integrable systems 
methods ) this theory remains to be rigorously constructed as an EQFT . 

Let us now consider a relatively simpler example, that of constructing a massive super-renormalizable 
scalar field theory. This has been studied in d = 3, with [<fi] = = |. We get £ = £, g = 
L~ N g, fi = L~ 2N £l and g is taken to be small. £ is marginal whereas g, \i are relevant parame- 
ters and grow with the iterations. After N iterations they are brought upto g, jl together with 
remainders. This realizes the so called massive continuum c/> 4 theory in d = 3 and this has been 
mathematically controlled in the exact RG framework. This is proved by Brydges, Dimock and 
Hurd in [BDH 1995] and earlier by Benfatto, Cassandro, Gallavotti and others in [BCGNOPS 
1980]. 

3 The exact RG as a Continuous Semigroup 

The discrete semigroup defined in (|2.13l of the previous section has a natural continuous counterpart. 
Just take L to be a continuous parameter, L = e* : t > and write by abuse of notation T t ,St, T t 
instead of T e t etc. The continuous transformations Tf 

T t F = S t fi rt *F (3.1) 

give a semigroup 

T t T s = T t+S (3.2) 

of contractions on L 2 (dfic) with fj,c as invariant measure. One can show that Tt is strongly 
continuous and therefore has a generator which we will call C This is defined by 

CF = lim r ^—^-F (3.3) 

whenever this limit exists. This restricts F to a suitable subspace T>(C) C L 2 (d^c)- con_ 
tains for example polynomials in fields as well as twice differentiable bounded cylindrical func- 
tions. The generator C can be easily computed. To state it we need some definitions. Define 
(D n F)(<fi; fi, . . . , f n ) as the re-th tangent map at <f> along directions fi,...,f n . The functional 
Laplacian Ap is defined by 

A r F(0) = J d Mr (C) {D 2 F){<t>- C, C) (3.4) 
where V = u. Define an infinitesimal dilatation operator 

V<J)(x) = x ■ V<p(x) (3.5) 

and a vector field X 

XF = - [4>\ (DF) (0; 0) - (DF) (0; V<j>) (3.6) 
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Then an easy computaion gives 



T t is a semigroup with C as generator, 
the linear PDE 



C = ±Aj, + X 
Therefore T t = e tc . Let F t (4>) 



(3.7) 

= T t F(4>). Then F t satisfies 



-i = CF > < 3 - 8 > 

with initial condition Fq = F. This evolution equation assumes a more familiar form if we write 
F t = e~ Vt , Vt being known as the effective potential. We get 

^ = CV t -\(y t ) 4) -{V t ) <t> (3.9) 

where 

• (W))* = y ^r(C)((W t )(</»; C)) 2 (3.10) 

and Vo = V. This infinite dimensional non-linear PDE is a version of Wilson's flow equation. Note 
that the linear semigroup Tt acting on functions induces a semigroup TZt acting non-linearly on 
effective potentials giving a trajectory Vt = IZtVo- 

Equations like the above are notoriously difficult to control rigorously especially so for large 
times. However they may solved in formal perturbation theory when the initial Vo is small via 
the presence of small parameters. In particular they give rise easily to perturbative flow equations 
for coupling constants. They can be obtained to any order but then there is the remainder. It 
is hard to control the remainder from the flow equation for effective potentials in bosonic field 
theories. They require other methods based on the discrete RG . Nevertheless these approximate 
perturbative flows are very useful for getting a prelimnary view of the flow. Moreover their discrete 
versions figure as an imput in further non-perturbative analysis. 



Perturbative flow : It is instructive to see this in second order perturbation theory. We will 
simplify by working in infinite volume ( no infrared divergences can arise because T(x — y) is of fast 
decrease ). Now suppose that we are in standard 4> 4 theory with [0] = and d > 2. We want to 
show that 

V t = jdx : |V^(x)| 2 : +g t : <P(x) 4 : : 0(x) 2 :) (3.11) 

satisfies the flow equation in second order modulo irrelevant terms provided the parameters flow 
correctly. We will ignore field independent terms. The Wick ordering is with respect to the 
covariance C of the invariant measure. The reader will notice that we have ignored a <p® term 
which is actually relevant in d = 3 for the above choice of [(f)]. This is because we will only discuss 
the d = 3 case for the model discussed at the end of this section and for this case the <fi 6 term 
is irrelevant. We will assume that are of 0(g 2 ). Plug in the above in the flow equation. 

\™' m : P nt m '■ represent one of the terms above with m fields and n derivatives. Because C is the 
generator of the semigroup T t we have 

o i\n,m 

( _ C )%' m : P n>m := (-^— - (d - mM - n)A"' m ) : P n , m : (3.12) 
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Next turn to the non-linear term in the flow equation and insert the (j) 4 term ( the others are already 
of 0(g 2 ) ). This produces a double integral of T(x — y) : cft(x) 3 :: </>(y) 3 : which after complete Wick 
ordering gives 

_^1 6 J dxdyt(x-y){: ^(xf^yf : +9C(x-y) : 0(x) 2 0(y) 2 : +3QC(x-y) : <j>{x)<j>{y) : +6C(x-y) 2 ) 

(3-13) 

Consider the non-local </> 4 term. We can localize it by writing 

: m 2 m 2 ■■= \ ■■ (m A + <p{yf - {m 2 - <p(y) 2 ) 2 ) ■ (3.14) 

The local part gives a cp 4 contribution and the the last term above gives rise to an irrelevant 
contribution because it produces additional derivatives. The coefficients are well defined because 
C, I 1 are smooth and T(x — y) is of fast decrease. Now the non-local 4> 2 term is similarly localized. 
It gives a relevant local (j) 2 contribution as well as a marginal |V0| 2 contribution. Finally the same 
principle applies to the non-local c/> 6 contribution and gives rise to further irrelevant terms. Then it 
is easy to see by matching that the flow equation is satisfied in second order up to irrelevant terms 
( these would have to be compensated by adding additional terms in V t ) provided 



(4-d)g t -ag 2 + 0(gl 
2m ~ hg 2 t + 0{gl 
cg 2 +0{gl) (3.15) 



dgt __ , 1 , 7w , _2 , r)ti ;>,, 
dt 

dt 

dt 

where a, b, c are positive constants. We see from the above formulae that up to second order in g 2 
as t — > 00, gt — > for d > 4. In fact for d > 5 the decay rate is 0(e - *) and for d = 4 the rate is 
0(t~ l ). However to see if Vt converges we have to also discuss the jiit,& flows. It is clear that in 
general the fit flow will diverge. This is fixed by choosing the initial fiQ to be the bare critical mass. 
This is obtained by integrating upto time t and then expressing fiQ as a function of the entire g 
trajectory up to time t. Assume that fit is uniformly bounded and take t — > oo. This gives the 
critical mass as 

f-OO 

fiQ = b I ds e~ 2s g 2 = fi c (go) (3.16) 
Jo 

This integral converges for all cases discussed above. With this choice of fio we get 



roc 



roc 

/ dse- 2s g 2 +t (3.17) 
j o 



'0 

and this exists for all t and converges as t — * oo. Now consider the perturbative £ flow. It is easy 
to see from the above that for d > 4, £t converges as t — > oo. 

We have not discussed the d = 3 case because the perturbative g fixed point is of O(l). But 
suppose we take in the d = 3 case [4>] = with e > held small as in [BMS 2003]. Then the 
above perturbative flow equations are easily modified ( by taking account of (|3.12|) ) and we get 
to second order an attractive fixed point g* = 0(e) of the g flow. The critical bare mass fio can 
be determined as before and the £j flow converges. The qualitative picture obtained above has a 
rigorous justification. 
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4 Rigorous RG analysis 



We will give a brief introduction to rigorous RG analysis in the discrete set up of Section 2 concen- 
trating on the principal problems encountered and how one attempts to solve them. The approach 
we will outline is borrowed from Brydges, Mitter and Scoppola [BMS 2003]. It is a simplification of 
the methods initiated by Brydges and Yau in [BY 1990] and developed further by Brydges, Dimock 
and Hurd [BDH 1998] . The reader will find in the selected references other approaches to rigorous 
RG methods such as those of Balaban, Gawedzki and Kupiainen, Gallavotti and others. We will 
take as a concrete example the scalar field model introduced earlier. 

At the core of the analysis is the choice of good coordinates for the partition function density 
z of Section 2. This is provided by a polymer representation ( defined below) which parametrizes z 
by a couple (V, K) where V is a local potential and K is a set function depending also on the fields. 
Then the RG transformation Tl maps (V,K) to a new (V,K). (V, K) remain good coordinates as 
the volume — > oo whereas z(volume) diverges. There exist norms which are suited to the fixed point 
analysis of (V,K) to new (V,K). Now comes the important point: z does not uniquely specify the 
representation (V, K). Therefore we can take advantage of this non- uniqueness to keep K small in 
norm and let most of the action of Tl reside in V. This process is called extraction in [BMS 2003]. 
It makes sure that K is an irrelevant term whereas the local flow of V gives rise to discrete flow 
equations in coupling constants. We will not discuss extraction any further. In the following we 
introduce the polymer representation and explain how the RG transformation acts on it. 

To proceed we introduce first a simplification in the set up of Section 2. Recall that the function 
u introduced before (|2.3|) was smooth, positive definite, and of rapid decrease. We will simplify 
further by imposing the stronger property that it is actually of finite range : u{x) = for \x\ > 1. 
We say that u is of finite range 1. It is easy to construct such functions. For example if g is any 
smooth function of finite range \ then u = g * g is a smooth positive definite function of finite 
range 1. This implies that the fluctuation covariance Tl of (|2.7jl has finite range L. As a result 
T n in (|2.10|) has finite range L n+l and the corresponding fluctuation fields Cn( x ) and £ n (y) are 
independent when \x — y\ > L n+1 . 

Polymer representation 

Pave R rf with closed cubes of side length 1 called 1- blocks or unit blocks denoted by A, and 
suppose that A is a large cube consisting of unit blocks. A connected polymer X C A is a closed 
connected subset of these unit blocks. A polymer activity K(X, 4>) is a map X, (j) — ► R where the 
fields (j) depend only on the points of X. We will set K(X, <fi) = if X is not connected. A generic 
form of the partition function density z(A, (j)) after a certain number of RG iterations is 

oo N 

-( A ) = E^ E e-^UK( Xj ) (4.1) 

N=0 ' X 1: ...,X N j=l 

Here the Xj C A are disjoint polymers, X = LiXj, and X c = A \ X. V is a local potential of the 
form (|2.23j) with parameters £, g, \x. We have suppressed the cj> dependence. Initially the activities 
Kj = but they will arise under RG iterations and the form ()4.1j) remains stable as we will see. 
The partition function density is thus parametrised as a couple (V, K). 
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Norms for polymer activities : 

Polymer activities K(X,(j>) are endowed with a norm ||i^(X)|| which must satisfy two properties: 

xnt = o^\\k 1 (x)k 2 (y)\\ < po || ||^i (y)|| 

\\T L K(X)\\<c^\\K(X)\\ (4.2) 

where X is the interior of X and \X\ is the number of blocks in X. c is a constant of 0(1). The 
norm measures ( Frechet ) differentiability properties of the activity K(X, 4>) with respect to the 
field <j) as well as its admissible growth in eft. The growth is admissible if it is jxc integrable. The 
second property above ensures the stability of the norm under RG iteration . For a fixed polymer 
X the norm is such that it gives rise to a Banach space of activities K(X). The final norm || • ||_4 
incorporates the previous one and washes out the set dependence 

||X|U = sup^ A(X)\\K(X)\\ (4.3) 

where A(X) = 7^( rf + 2 )l^l . This norm essentially ensures that large polymers have small activities. 
The details of the above norms can be found in [BMS 2003]. 

The RG operation map / is a composition of two maps. The RG iteration map z — > T^z in- 
duces a map V — > V L and a nonlinear map Tl ■ K — » K = Tl(K). We then compose this with a ( 
nonlinear ) extraction map 8 which takes out the expanding ( relevant) parts of K — > S(K) = K' 
and compensates the local potential VI — > V such that T^z remains invariant. We denote by / 
the composition of these two maps with 

V^V' = f v (V, K), K^K' = f K (V, K) (4.4) 

The map Tl : 

Consider applying the RG map Tl to l|4.1j) . The map consists of a convolution /Ur L * followed 
by the rescaling Sl- In doing the integration over the fluctuation field Q we will exploit the inde- 
pendence of C(x) and ((y) when \x — y\ > L. To do this we pave A by closed blocks of side L called 
L-blocks so that each L-block is a union of 1-blocks. Let X L be the L-closure of a set X, namely the 
smallest union of L-blocks containing X. The polymers will be combined into L-polymers which 
are by definition connected unions of L blocks. The combination is performed in such a way that 
the new polymers are associated to independent functionals of £. 

Let V(X,cp)), to be chosen later, be a local potential independent of £. For the coupling con- 
stant sufficiently small there is a bound 

He-^ll <2l y l (4.5) 
We assume that V is so chosen that the same bound holds when V is replaced by V. Define 

P(A, C, 0) = e~ v ^ + ^ - e -^( A <<» (4.6) 
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Then we have 

e -v{x c ,^) = e -v(x c , f +0) = JJ {e -v(&,4>) + p(Aj c 0)) (47) 

Aclc 

where X c is the closure of X c . Expand out the product and insert into the representation (|4.1|) for 
z(A, £ + </>)). We then rewrite the resulting sum in terms of L-polymers. The sum splits into a sum 
over connected components. Define for every connected L-polymer Y 

1 _ N M 

BK{ y)= E ww S ^ (Xo) n^i)II^) ( 4 - 8 ) 

AT+M>1 ' ' (X,-),(Ai)-*Y 3=1 i=l 

where Xo = Y \ (UXj) U (UAj) and the sum over the distinct Aj, and disjoint l-polymers Xj is 
such that their L-closure is Y. (|4,lj) now becomes 

-( A ) = E^ E e-^n^O ( 49 ) 

' Yi,..,Yjv 3=1 

where the sum is over disjoint, connected closed L-polymers. We now perform the fluctuation 
integration over £ followed by the rescaling. Now V(Y C ) is independent of £. The £ integration sails 
through and then factorizes because the Yj being disjoint closed L-polymers they are separated 
from each other by a distance > L. The rescaling brings us back to l-polymers and reduces the 
volume from A to L _1 A. Therefore 

1 - N 

z'(L- 1 A) = T L z(A)=^2 m Yl e- v ^ H(T L BK)(X,) (4.10) 

' Xi,..,X N 3=1 

where the sum is over disjoint l-polymers, X c = L _1 A \ X. By definition Vl(A) = SlV(LA) and 
(TlBK)(Z) = Sl^Fl * BK(LZ). This shows that the representation (|4.1[) is stable under iteration 
and furthermore gives us the map 

V ^ V L 

K^K = f L (K) = T L BK (4.11) 

The norm boundedness of K implies that Tl(K) is norm bounded. We see from the above that 
a variation in the choice of V is reflected in the corresponding variation of K. The extraction 
map £ now takes out from K the expanding parts and then compensates it by a change of Vl 
in such a way that the representation Q4.1UJI is left invariant by the simultaneous replacement 
Vl — > V, K — > K' = £{K). The extraction map is non-linear. Its linearization is a subtraction 
operation and this dominates in norm the non-linearities, [BDH 1998]. 

The map V — * Vl — * V' leads to a discrete flow of the coupling constants in V. It is conve- 
nient to write K = K pert + R where R is the remainder. Then the coupling constant flow is a 
discrete version of the continuous flows encountered in Section 3 together with remainders which 
are controlled by the size of R. In addition we have the flow of K. The discrete flow of the pair 
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(coupling constants, K) can be studied in a Banach space norm. Once one proves that the non- 
linear parts satisfy a Lipshitz property the discrete flow can be analyzed by the methods of stable 
manifold theory of dynamical systems in a Banach space context. The reader is referred to the 
article [BMS 2003] for details of the extraction map and the application of stable manifold theory 
in the construction of a global RG trajectory. 

5 Further topics 

1. Lattice RG methods: Statistical mechanical systems are often defined on a lattice. Moreover the 
lattice provides an ultraviolet cutoff for euclidean field theory compatible with Osterwalder-Schrader 
positivity. The standard lattice RG is based on Kadanoff- Wilson block spins. Its mathematical 
theory and applications have been developed by Balaban starting with [B 1982] and Gawedski and 
Kupiainen starting with [GK 1980]. This leads to multiscale decompositions of the Gaussian lattice 
field as a sum of independent fluctuation fields on increasing length scales. Brydges, Guadagni and 
Mitter [BGM 2004] have shown that standard Gaussian lattice fields have multiscale decomposi- 
tions as a sum of independent fluctuation fields with the finite range property introduced in section 
4. This permits the development of rigorous lattice RG theory in the spirit of the continuum frame- 
work of the previous section. 

2. Fermionic field theories: Field theories of interacting fermions are often simpler to handle 
than bosonic field theories. Because of statistics fermion fields are bounded and perturbation series 
converges in finite volume in the presence of an ultraviolet cutoff. The notion of studying the 
RG flow at the level of effective potentials makes sense. At any given scale there is always an 
ultraviolet cutoff and the fluctuation covariance being of fast decrease provides an infrared cutoff. 
This is illustrated by the work of Gawedzki and Kupiainen [GK 1985] who gave a non-perturbative 
construction in the weak effective coupling regime of the RG trajectory for the Gross-Neveu model 
in two dimensions. This is an example of a model with an unstable Gaussian fixed point where the 
initial coupling has to be adjusted as a function of the ultraviolet cutoff consistent with ultraviolet 
asymptotic freedom. 
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